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Abstract - Coupling Inductors between two channels can 
improve both rat steady state and dynamic performances of 
ranlti-caaanei tatarteaviog >oitagc regulator ooodiies (VRMs). 
However, the complicated electrical circuit node! for the muhi- 
eoopling core structures prevents the coupling concept from 
being extended to integrated muirJ-channeJ core structures. This 
paper presents a novel circuit modeling concept that can be 
•asUy applied to different complicated multi-coupling core 
structures. Compared to conventional modeling methods, the 
proposed method involves fewer Inductors. It is more directly 
related to the physical core structures and Is easier to converge 
in computer simulation. Simulation, results are presented zo 
verify the validity of this modeling coacepi. 

L Introduction 

Multi-channel interleaving voltage regulator modules 
(VRMs) have been widely adopted as industry practice. With 
proper coupling design, the integration ox" the inductors in cwo 
of the interleaving channels with 1 SO 3 phase shift inrproves 
both the steady-state performances and the transient 
responses [1]. The coupling effects reduce the steady-state 
current ripple and :he equivalent transient inductance. Tne 
conduction losses in the MOSFETs can thus be reduced. The 
integrated coupling core structure also makes the fluxes more 
evenly distributed in the core, winch reduces core losses. 
Experimental results show chat the efficiency irrrprovement 
due to the coupling core structure is significant. Moreover, 
the integrated coupling core structures rnaks the manufacture 
of the cores easier and improve their mechanical stability. 

For increasing current ratings of the VRMs, more 
interleaving channels become more desirable. Current 
research interest focuses on the use of integrated coupling 
core structures to improve the performance and increase the 
power density of multi-channel applications. However, the 
defects cf existing electrical circuit modeling methods for 
complicated cere structures prevent various studies of the 
different integrated coupling core structures. These defects 
induce complexity, involving many inductors and controlled 
sources, converging problems in computer simulations, 
elements without clear physical meaning, and lack of 
scaiabilities. This paper proposes a solution to these problems 
by presenting a novel electrical circuit modeling concept for 



the magnetic core structures. The proposed modeling concept 
uses minimum inductors and is easier to converge. The 
proposed modehng concept is closely related to the magnetic 
reluctance model and the physical structure of the ceres. The 
proposed modeling concept is easy to scale to different 
coupling channels. 

H. Review cf Existing Methods 

The modeling of the integrated muln-chanr.cl core 
structures is similar tc chat of the muiti- winding rransfenners. 
unless ±e leakage inductance in the transformers is much 
smaller than the magnetizing inductance. The leakage 
inductance in the transformers is actually the self- inductance 
in coupling inductor structures, while the common 
magneriang inductance in transformers is the coupling 
inductance in coupling inductor structures. Modeling of the 
multi-winding transformer is an extensively researched area. 
This section reviews the sxisting common methods. 

For two-winding transformers, the 7>model is the most 
commonly used The r:-mcdel car. he easily simplified and 
applied in rwo-channei coupling inductor sttuctutcs [1], as 
shewn in Fig. I. One extension of the z-model fcr three- 
channel coupling inductor structures is shown in Fig. 2. 
However, this model is lot a general method. It is only valid 
for the special case in which the couplings among the 
channels are ail the same. This cannot always be the case. As 
will be discussed later, some s> f nnnetric core structures may 
also result in different couplings among channels. Tne 
general madiernarical model for multi-channel coupling 
inductor cases is described as Formula 1* There is a total cf 
n(n-H)/2 independent components in the n^n diagonal 
symmetric ma irix. In order tc describe the case generally, the 
circuit model needs to include n(n-iV2 inductors in an n- 
channel coupling structure. A more general extension of the 
rr-rnodei for three -channel coupling inductor structures is 
given in Fig. 3 [3). which shows that this model is actually 
difficult to ippiy to coupling Indnrjor structures because of 
the different defiriuon of terminations and parameters in the 
model. Moreover, this model is also difficult tc extend to 
more than ±ree channels. 
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Fig. 



. Simplified n-rnodd for two-channel coupling inductors. 




Fig. 2. A special extension of r-mocel for three-channel coupling inductors. 




Fig. 3. A general extension of z-nwdei for Lhrec-.^inncI ccupJmg inductee. 



(i) 



Considering the model shown in Fig. 1, the star connection 
of the three inductors can be replaced by the delta connection 
of the other three inductors with equivalent transformation, as 
shown in Fig. 4. The delta connection can be extended to 
multi-channel coupling structures. The delta connection 
modeling of a four-channel coupling structure is shown in 
Fig. 5 [4], However, the inductor? in the model are 
mathematical transformation results without clear physical 
meanings. The derivation and explanation of the parameters 
in the model are very difficult 
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Fig. 4. Y-A trajisforraation of x-mode; for Two-channel couplirg inducttn. 
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Fig. 5. Delta ccnnccron model for ihree-ctamet couplir.8 inductors. 

The general procedure to derive the electrical circuit model 
starts from the magnetic reluctance model. An example of a 
four-channel coupling inductor structure and its reluctance 
model are shown in Fig. 6. The magnetomotive force (MMF) 
sources, which are the products of *iie currents and turns of 
the windings, translate into voltage sources in the reluctance 
model to represent the windings, as shewn :n Fig. 5. The 
reluctance model, which can be easily derived from the core 
geometry and the material permeability, directly describes the 
core. The paramr.srs x the mocei have clear physical 
meamngs. However, this mod*! cannot bs directly used in 
circuit simulation to ascertain the electrical effect of the 
magnetic components. The reluctance mode: of the core 
needs to be transferred to the electrical model which is the 
difficult pan or the model mg. The major ^xisrutg methods are 
as follows. 




Fig. 5. Derivation of r5iucjr.ee rrcicl f-o.T :ore 5^;-rjre. 




Fig. 7. A general irodiiicaior of >noce: :> -.hrse^ar.r.el coupiing ;cr= 
itructire. 

One method modifies the model shown m Fig. 2. The 
modification is shown in Fig. 7 [5]. The couplings between 
the windings are represented ii the current sources and the 
inductance. A total often inductors 2nd six centre lied current 
sources ere used for tarce-whaziurl coupling inductors. This 
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model is apparently :oo complicated and ;a converging 
problems in computer siniaiarioa maxe :i jven more 
pro b Icinsnc . 

Another method includes me magnetic reluctance mocei 
directly is the electric circuic simulation [6]. An mrffrfacrr 
between the magnetic and the eiecrric circuit is accessary to 
transfer the parameters between the magnetic circuits and the 
electric circuits, as shown in Fig. 3. The interface relates the 
current and voltage in the electric circuit with me flux and 
MMF source in the magnetic reluctance model. Although the 
concept of the interface is simple, the implementation 
involves cnnent-controilec current sources, current- 
controlled voltage sources* and inductors. The complexity of 
the interface results in rime -consuming simulation and 
converging problems. 




Fig. 3. Including r?!«ctar.oc n>^ltl ;n circuit simulation 



The last method applies duality rjansformaticn rules to 
generate the electric model from the reluctance model 6]. 
However, for complex core structures, the transformation is 
labor-iniensive and has high probability of error. Smce the 
model requires ideal transformers, it also has converging 
problems in circuit simulation. 

The review of the existing modeling methods shews that 
they all have defect, such, as cc-nmiexiry. converging 
problems in simulation, and difficulties in racdel derivation 
and explanation. In order to fine a proper modeling method 
for multi-channel coupling inductor structures, this paper 
proposes a novel racdeHng method that snns from a di&reni 
reluctance modeling concept. 

m. Concspto? The Proposed Method 

in Che magnetic reluctance model the reluctances are only 
determined by core structure and material permeability. I: is 
straightforward and cannot changed. The windings are 
represented as the MMF sources, which art the currents that 
run through the windings. The existing modeling methods all 
use this method to derive the reluctance model 

From another point cf view. the voltages applied to the 
windings are also die parameters that can represent the 
windings. According :c Fatacay's Law. the vclnge applied to 
the winding integrated with the time {whici ts mere simply 
stated as the voitage-secend of tie wxdicz} determines the 
flux in the corresponding core '.eg, as shown m Formulas 2 



and 2. The dux in the core is represented as the current in the 
reluctance model. Thus, the windings can also be represented 
as the current sources m the reluctance model. The values of 
the current sources are determined by the voltage* second of 
the windings. To distinguish this reluctance model rxom the 
conventional reluctance model, this one is named the flux 
source reluctance model while the conventional model is the 
MM? source reluctance model. The two models are 
equivalent ?cr the MMF source model the MMFs are the 
known parameters, and the fluxes in the branches need tc be 
determi ned. For the flux source model, the fluxes in die 
branches are the known parameters, and the MMFs of the 
loops aeed to be determined Figure 9 is the flux source 
reluctance model for the core structure shown in Fig. 6. The 
only differences berwesa the two reluctance models ire the 
sources that represent the windings. 

N-0~\y>dt (2) 



dt 
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Rg. 0. Fi'jx source reluctance *vcaz- :cr 6c care str^ct'jre f !«. S 

For the flux source 'reluctance model, the MMF of a certain 
loop is cetcrmmed by its relutances and the fluxes through 
them. For the flux source reluctance model, the fluxes m the 
branches become very clear. The values can be easily gotten 
solely from observation Accordingly, the MMF of each loop 
can also be determined easily. However, for the MMF source 
model, it is usually necessary to solve a set of a coupled 
equations m order to 6nd out the branch fluxes in an n- 
channel coupling cere structure. For the modeling rrte ±ods 
shown in Fig. 2 and Fig. 5. so King the equation set is done by 
hand, which mzJ&es the parameters in the model difficult to 
derive and without clear physical meaning. For the other 
modeling methods reviewed in Section U. the equation set is 
solved in the circuit simuiarion. which makes the simulation 
rirrc-censumirig and causes convergence problems. In other 
words, the problems cf the existing modeling methods arc 
caused by the difficulty of solving the coupled equations. 

The flux source reluctance model also requires solving a 
set of c equations, but the equations are decoupled. Solving n 
decoupled equations is much simpler than solving n coupled 
equaacrs, especially when n increases. Although me flux 
source reluctance mode: ind the MM? source reluctance 
model are equivalent, solving the flux source reluctance 
model is much suqrlcr This is the main reason that the 
circuit model based on the ilux source reluctance model is 
much simpler than the circuit mode! based on ±e MMF 
source reluctance model 

For the model shewn in Fig. 9> a general cell can be 
extracted, as shown in Fig. 10. According to Airperc's Law. 



101 



the MMF of Loop \ which represents the cmrern; in the 
winding, can be described in formula 4. The forma la can be 
very easily derived jusi by observing Fig. 1 0. 




Fig. 10. A general celt in the flux source reiucance node]. 

The sum of the reluctances in Loop i, 
Rt =S^J - ^^t)' * defined as the self- 

j 

reluctance of Loop L Tie reluctance involving both Loop i 
and j, % , is defined as the rnnruai reluctance between these 

two loops. As shown in Forrmia rhe MMF can be divided 

into two pans: the self-MMF, which is generated by flux *j> { 

through the sdf -reluctance of the loop, and ".he mutual MMF, 

which is generated by fluxes otasr than 4> ; through the 

corresponding mutual-rcfcctances. 

Formula 5 is the differentiation of Formula 4. Submitting 
Formula 3 into Formula 5 results in Formula 6. Formula 6 
represents the current through the corresponding winding. 



di 



dt 



as 



(5) 
(6) 



The definition of the inductances is shown in Formula 7. 
The inductance L u is denned as the seU-bductance of Loop 
i, which corresponds :o ±e self-reluctance of ±e loop. The 
inductance L UJ is defined as the mutual inductance between 

Loop i and j, which corresponds to the mutual reluctance 
between the two loops. All these inductances have clear 
physical meanings, and are closely related to the magnetic 
core structure. 



j 

V -JV 



self -inductance 



mutual -inductance 



(7) 



With the definitions in Formula 7, Formula 6 can be 
simplified into Formula S. 



**** _ 
<ir L 



1 



(S) 

The current through a winding can be divided into two 
parts; the 3cif-currenl> which is the firs: item in Formula 8, 



and the mutual (cr coupling) currents, which is the rest of the 
items in Formula 8. The definitions of self-current and rnutuai 
current are shown in Formula 9. Formula 3 can be rewritten 
as Formulas 10 and 11. 
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mutual-current 
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The voltage applied to die winding and die self- inductance 
determine the self-current through the winding. The voltage 
on the other windings and the corresponding mutual 
inductances determine the mutual currents. The self-current is 
die decoupled pan of :he winding current. The mutual 
currents are the coupling effects among different windings. In 
this specific case, the mutual currents have only two parts, 
which correspond to the two adjacent windings. In physical 
reality, there are coupling effects among all the windings. 
This model has simpler expression but the coupling effects 
remain. As it is shown in Section IV, this model does result in 
the coupling of all the windings, as da the ether models. 

The total current through a winding is generated by several 
winding voltages. In other words, the voltage zpphed to a 
certain winding not oniy generates the sclf-currcn: a the 
winding, bur also generates the mutual currents in the other 
windings. The mutual and self-currents generated by the 
same voltage have waveforms proportional to each ether. The 
mutual-current waveforms can be described as the self- 
current waveform generated by the same voltage with a 
proportional coefrkient. Another format cf Formula 9 is 
shown in Formula 12. The coefficient between the mutual 
and self-currents generated by the same voltage can be found 
by comparing the two equariens in Formula 12, The result is 
shown in Formula 13. The coefficients a {J are oriiy 

determined by the reluctances. As previously discussed, the 
parameters are directly rented to the core structure, . can be 
easily derived, and have clear physical meanings. 

_ Iv* 



t. , = 
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(13) 



With Formula '3, Formula *. I can be rewritten as FcmuLa 



14. 
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h=hx' 



The mamx format of Formula 14 is shown in Formula 15. 
The total currents in the windings -an be represented as the 
linear transformation of the seif-cuirems in cite windings. 
Matrix A in Formula 15 represents the linear transforroatioa 



I 



j L'«J 



Ha 
hi 



As discussed previously, the self-currenr in a winding is 
the decoupled pan of the wincing anient. It is determined 
only by the corresponding individual winding voltage. 
Formula 15 can be further written as Formula 16. Formula 16 
is a general expression of the results of the flux source 
reluctance model for n-chanad coupling magnetic structures. 
For known voltage sources applied to the corresponding 
windings, Formula 16 gives the current waveform in the 
winding The required parameters in Fonnuia 16 are the 
coefficients of matrix A and the self-inductances L tJ . As 

discussed, all the parameters in the formula can be derived 
directly from the reluctance model solely by observation. The 
results can be reached by a simple maocix multiplication 
avoiding problematic solving of n-coupled equations 
encountered in conventional modeling methods. Thus, this 
method is much simpler. 



(16) 




From another point of view, Formula 1 is the general 
expression of the results of the MMF source reluctance 
model. Formula 6 is just the inverse transformation of 
Formula 1. Matrix A in Formula 16 is the inversive mamx of 
the inductance matrix in Formula 1. The elements of Matrix 
A in Formula 16 are cirec-Jy related to the reluctance model 
with clear physical meanings, while the elements of the 
inductance matrix in Formula 1 cornet be easily related to the 
reluctance model To derive the parameters in Formula 1, 
matrix inversion is unavoidable Matrix inversion is actually 
the solving of the linear equation set mentioned previously. 
The current source reluctance model avoids matrix inversion, 
thus reducing the complexity. 

Toe next step is to Implement the model in circuit 
simulation. From Formula 16. the model requires n 



independent inductors for n-channci coupied core structures. 
The elements in Matrix A. a : ; . , can be represented by the 

controlled current sources. The simu l ation model for ihe 
general ceil shown in Fig. 10 and Formula 6 is shown in Fig. 
It. The ftmcaxm of die controlled current sources in Fig. 1 1 
is to add up the currents in the different inductors. They arc 
not involved in circuit simulation as (bey are in the model of 
Fig. 7. Thus, the proposed method does not have converging 
problem seen in the conventional methods' simulations. 




Fig. 1. Sirttlation mocd for :hc gersenl ceil. 

iv*. ver5icatign aind sumuut.on results of the 
Proposed Model 

In order to verify me validity of the proposed modeling 
method, simulation circuits with coupling cere structures are 
built based or. both die proposed model and the conventional 
models. The circuit is an n-channei interleaving quasi-sqaare- 
wave (QSW) convener with coupling core structures and is 
shewn in Fig. 12. As mentioned previously, for two-channel 
coupling core structures, the conventional magnetic modeling 
(shown in Fig. 1) is sufficiently accurate and fast. For two- 
channel coupling core structures, the proposed magnetic 
model will compare with the model shown in Fig, 1 . 
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F:g. 12. Imritaving QS'»V soirHtfter "dth coupling core itrucruttS tor 

The simulation rcsuks of the winding current for both 
models are shown in Fig. 13 by replacmg the shaded area in 
Fig. 12 with the model shown in Fig. 1 and It. The two 
modeling methods give exactly the same results. 
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(a) conventional model 
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(b) proposed modef 
rig. 13. Siiruiition results for two-channel hterleavrrg QSW converter with 
co-aping core structures. 




5ig. 14. Reluctance mods! of * two-charmel cnupling core structures. 



through the sdf-inencters in Fig. 11. The controlled current 
sources in the proposed model (as shown in Fig. 1 1) are just 
to get the weighted sum of the flux waveform. The winding 
current can be easily gotten using the proposed method. Thus, 
there is no converging problem in thecircuit simulation. 

Hie channel current simulation results for three-channel 
coupling core structures for the conventional model and the 
proposed model are compared in Fig. 16. The two simulation 
results are exactly the same. These simulation results confirm 
the validity of the proposed model for coupling core 
structures. 



i a) conventional modei 




Fig. 15. Analysis of the waveforms in the proposed method. 

For the integrated magnetic core structure shown in [;], the 
reluctance model of the core is shown in Fig. 14. As 
illustrated in Formula 6, the proposed magnetic model is 
actually a weighted sum of the different voltage-seconds 
applied to die integrated magnetic core. To explain the 
concept more clearly, the process is illustrated in Fig. 15. For 
the two-channel interleaving QSW converters, the voltage 
waveforms applied to the two windings arc the waveforms of 
vl and v2 shown in Fig. 15. The integranon of me waveforms 
gives the flux waveforms in the corresponding core legs. The 
flux waveforms wi:h the weighted coctrtecnt arc the 
waveforms of (R^R f )^ x and R f • £ in Fig. 15. The 
current waveform in the winding, which is the waveform of 
il in Fig. 15, is just die sum of the flux waveforms of 6, and 

p 2 with weighted coefficients of (R y + R e ; and R : , 
respectively, Toe flux waveforms are actually proportional to 
the self-current waveforms, which are also the currents 



(b)proccsed mccel 
Fig. \6. Sirmilarion results fcr three-channel Interleaving QSW :*nv«*cr 
w.th coupling core 3trucrures. 

As mentioned previously, the purpose of mis model is to 
be able to be applied to more complicated coupling core 
structures in order to identify the performances of different 
core coupling structures. Figure 17 gives the simulation 
current waveforms of the four-chancel integrated coupling 
core structure shewn in Fig. 6 with a certain winding 
arrangement. The detail pertbrmar.ee evaluancns cf the 
different coupling core structures are beyond the scope of this 
paper. The purpose of this paper is to show mat the proposed 
magnetic model can easily handle a variety of complicated 
coupling core structures. 



Fig. 17. Simulatior Jesuits for riw ceupltn$ cert structure in Fig. 6. 
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For the specific core structure shown in Fig. 6, Matrix A in 
Fonmila 16 fass the following fowidt 



V* V V o 

0 0 



-1 



X 4J ^44 J 



(17) 



From the viewpoint of Ac flux source modeling, the fluxes 
are only affected by the self and ±c two adjacent voltage 
sources. That is the reason for the zeros in the matrix Of 
course, in reality all the windings are coupled. The flux 
source modeling does not change the feet but simplifies the 
model. 



M 3 



(l 2 m) 



The inductance matrix for the conventional magnetic 
modeling in Formula I is the inversion of Matrix A, which 
could be overwhelmingly complicated For simplification, 
ayniTre that aft the mutual inductances in Formula 17 are the 
same as M The inversion of Matrix A is shown in Formula 
18. It is much more complicated than Formula 17. This 
implies that the proposed flux source magnetic modeling 
method gives much simpler simulation models for multi- 
channel coupling core structures. From Formula 18, it is also 
clear that all the four inductors are coupled Although the 
coupling is not explicidy shown In Formula 16, die two 
formulas describe the same fact. The simulation waveforms 
and the formulas also show that the seemingly symmetric 
core structure actually gives an asymmetric coupling pattern. 
Only with the accuracy and simplicity of the proposed 
magnetic modeling concept can the analysis of these subtle 
performance differences become possible. 



-1 vM'^LV 1 \M--L-j -1 
» L 3 ' M 1 
•1 



1 

(m»-0 



1 vIVT-L*' -1 \M'-L1 



M 2 



I -t (m 2 -l 2 ) [\r-2-L 2 ) 



M 3 



(IS) 



V. Conclusion 

The conventional magnetic modeling methods are too 
complicated to be useful for multi-channel coupling core 
structures. This paper presents a novel flux source modeling 
concept for the magnetic modeling of multi-channel coupling 
core structures. The method is much simpler than the 
conventional methods. Moreover, the proposed method is 
directly related to the magnetic reluctance model Every 
elemenr in the model has. clear physical meanings. Circuit 
simulation results for both methods are compared to verify 
the validity of the proposed method An example of a four- 
channel coupling core structure is provided to show the great 
reduction of complexity achieved by the proposed modeling 
method compared to the conventional methods. 

The proposed magnetic modeling method is suitable to 
evaluate the performances of die different multi-channel 
coupling core structures. 
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